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We theoretically consider specific adhesion of a fluctuating membrane to a hard substrate via the 
formation of bonds between receptors attached to the substrate and ligands in the membrane. By 
integrating out the degrees of freedom of the membrane shape, we show that in the biologically rel- 
evant limit specific adhesion is well described by a lattice gas model, where lattice sites correspond 
to bond sites. We derive an explicit expression for the effective bond interactions induced by the 
thermal undulations of the membrane. Furthermore, we compare kinetic Monte Carlo simulations 
for our lattice gas model with full dynamic simulations that take into account both the shape fluctu- 
ations of the membrane and reactions between receptors and ligands at bond sites. We demonstrate 
that an appropriate mapping of the height dependent binding and unbinding rates in the full scheme 
to rates in the lattice gas model leads to good agreement. 

PACS numbers: 



I. INTRODUCTION 

Adhesion of biomembranes to each other is ubiquitous 
in living organisms since it is inter alia crucial for pro- 
cesses like cell signaling or wound healing [T]. While in 
many cases two membranes may effectively be attracted 
to each other due to non-specific interactions, like electro- 
static or van der Waals interactions, adhesion in biolog- 
ical systems is typically supported by receptors and lig- 
ands in the membrane that form bonds upon contact [2]- 
[5]. The largest protein families involved in adhesion are 
cadherins, integrins, and selectins [6 . In order to gain 
understanding of the processes of specific adhesion much 
experimental effort has been made to develop model sys- 
tems using artificially prepared lipid bilayer vesicles with 
inserted ligands that are brought into the vicinity of re- 
ceptors tethered to stiff or soft, polymer-cushioned, sub- 
strates [THTT]. 

From a theoretical perspective, membranes are well de- 
scribed as two-dimensional sheets with a bending rigidity 
and a rather small effective surface tension. This model, 
first introduced by Helfrich [T2| , has been successful in de- 
scribing morphology and dynamics of free lipid bilayers 
and vesicles |13| . Since membranes are often not free but 
move close -or adhere- to other membranes (or, in model 
systems, to substrates) much work has been dedicated to 
the understanding of non-specific interactions between 
membranes |14| . In particular, analogies to wetting the- 
ory and the character of the transition between bound 
and unbound membrane [15H18j have been explored. In 
many studies on adhesion receptor-ligand bonds were in- 
tegrated out exactly in order to derive an effective inter- 
action potential that depends on the distance between 
the membranes or between membrane and substrate |19l - 
[55]. Other studies employed mean-field approaches that 
use coarse-grained bond density fields [551 [51] . In all of 
these approaches the discrete nature of specific adhesion 
involving receptor-ligand pairs is hidden. It only remains 
visible in more recent work relying on the use of simu- 



lations that typically combine the continuous nature of 
the membrane with the discrete nature of the adhesion 
bonds M.21,- 

In this work, we follow the opposite strategy com- 
pared to previous theoretical studies by integrating out 
the membrane shape in order to describe adhesion purely 
through receptor-ligand bonds. Starting from a mixed 
continuous-discrete model we show that in the regime of 
confined and rather stiff membranes, which is relevant for 
cells, it maps onto a lattice gas, where lattice sites corre- 
spond to receptor-ligand pairs. Membrane properties en- 
ter the effective interactions between bonds. The benefits 
from such an approach are two-fold: firstly, we can ap- 
ply the vast knowledge and methods developed for lattice 
gases to the adhesion behavior of membranes, and, sec- 
ondly, it allows for the efficient numerical simulation by 
only taking into account the discrete bonds. We demon- 
strate in which regime such a pure lattice gas description 
is valid and explicitly calculate the effective bond inter- 
actions. By constructing realistic reaction rates taking 
into account the distance between the membrane and the 
substrate, we are capable of qualitatively capturing the 
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FIG. 1: Sketch of the membrane. Besides non-speciflc inter- 
actions between substrate and membrane, bonds (with rest 
length lo) between receptors and ligands embedded into the 
membrane can form. The mean height of the membrane with 
respect to the substrate is h. Thermal fluctuations of the 
membrane are of size (jq. Above bond i, the deviation of the 
membrane height profile from the mean height is aohi. 
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dynamics of adhesion within the lattice gas framework. 
We support this claim by comparing simulation results 
for the full Hamiltonian with lattice gas simulations em- 
ploying the kinetic Monte Carlo method. 



II. MODEL AND MAPPING 

We consider a fluctuating bio-membrane in a fluid at 
temperature T above a substrate. The membrane is de- 
scribed in the Monge representation by the height profile 
/i(r) with r — {x,y). The projected area is A and we 
employ periodic boundary conditions. In addition, N re- 
ceptors at positions {r^} are embedded in the substrate 
that can bind to the corresponding ligands in the mem- 
brane, see Fig. [T] 

Our starting point is the total Hamiltonian H = Ho + 
T^ns+Hs composed of three terms [201 113 ■ The first term 
is the Helfrich energy 



no[h{r)] = 



dr^[V^Mr)]^ 



(1) 



with bending rigidity k governing the thermal fluctua- 
tions of the membrane. Non-specific interactions between 
the substrate and the membrane are due to, e.g., sur- 
face charges, van der Waals forces, or steric repulsion. 
They occur over ranges where individual molecules are 
not resolved. These non-specific interactions are mod- 
eled through the simple quadratic potential 



"Hi, 



[h{r)] = [ 

J A 



(2) 



with strength 7 resulting from a Taylor expansion around 
the minimum at height ho- The magnitude of the mem- 
brane height fluctuations in the absence of bonds is 



^ mO) - h] 



dq 



1 



(3) 
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using 



Eq. (Al) for the deflnition of the Fourier coefficients). 
Here, h is the mean height of the membrane. 

The focus of this work lies on specific substrate- 
membrane interactions through ligand-receptor pairs 
contributing with 



N 



[/i(r) ^ fej -(Ji + ao/i« - 



(4) 



to the total Hamiltonian. This energy is determined by 
two competing terms, the specific binding energy Cb and 
the energy stored in a stretched bond with effective spring 
constant k. The N binary variables bi are 1 for a closed 
bond and for an open bond. The bonds have rest length 
Iq, and h{ri) = h + (J^hi is the distance of the membrane 



ligand at Vi from the substrate. The dimensionless hi 
are of order unity. It will be advantageous to deflne two 
more quantities: the (average) area per bond 



and the fraction 



A/N 



1 ^ 



N 



1=1 



of closed bonds. 

The partition sum reads 



Mr)] 



-mih{r)-{b.}] 



(5) 



Our goal is to derive an effective lattice gas model by 
integrating out the height fleld h{r). The integration 
over the membrane fluctuations under the constraint of 
given heights hi is carried out in appendix \K\ 



E 

{b^} 



dhi ■ ■ ■ dhN e 



-hi:.,{rn-^).,h.h,~l3E^ 



(6) 



leading to a coupling of the height variables hi with cou- 
pling matrix 



cosq • (rj - Vj) 



(7) 



Kq^ + 7 

where the sum runs over independent wave vectors only. 



We rewrite the explicit expression for from Eq. ( A2 ) 



in terms of three dimensionless parameters Xj and A: 
l[xJ:^b^h^-N9X]^ 



NO + N(l)x 



1 



N 



-ex' 



i=l 



[(3A{Kq^ +7)]"^ with p EE l/ksT (see Here, 



(8) 



(9) 



is an effective cooperativity parameter between bonds 
discussed further below, while 



6 = p-ia^al 



(10) 



is the average energy contribution from the non-speciflc 
potential per bond. Finally, A = (/iq — lo)l<^a is the ratio 
of the distance between the minimum of the non-speciflc 
potential and the rest length of the tether, and the am- 
plitude of height fluctuations. 

For now we assume that x is small. Below and in the 
next section we will demonstrate that this assumption 
indeed corresponds to the range of parameters we are 
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interested in. Hence, Taylor-expanding Eq. ([8]) up to 
first order in x we obtain 



N 



N 



with the typical tether energy 



Xhl - /3(eb - et) 



(11) 



Within this approximation becomes independent of 9. 
Instead of A we will take et as the independent variable 
in the following, leaving us with et, Cb, and x- We now 
perform the Gaussian integration in Eq. ([6|. To this end, 
we have to invert the matrix {m'~^)ij +xbiSij, leading to 
rriij + 0{x)- Hence, to first order in x we finally arrive 
at 



(12) 



which is isomorphic to the lattice gas model of the liquid- 
gas transition with a binary density bi |28j . The sum runs 
over all bonds where the diagonal term is included in the 
lattice gas chemical potential 




= eb - et(l - x)- 



(13) 



The effective interaction energy between bonds is given 
through 



i^ij = etx^dr^ 
with length scale 



,1), m(r) 



e^(^/7)i/^ 



Ik ctq 



'-keio(r/0 



(14) 



(15) 



For the derivation and form of this effective interaction, 
see appendix |B] and Fig. [2] 

This mapping from the original full model to an ef- 
fective bond lattice gas is the central result of this paper 
valid for any geometry of bonds. The small quantity x de- 
fined in Eq. (|9| determines the local cooperative interac- 
tions of bonds. The physical picture is that a closed bond 
pulls down the membrane, assisting neighboring bonds to 
form. We have assumed the local effect of this behavior 
to be small; either through weak links (small k) such that 
the deforming force on the membrane is small, or by a 
stiff, confined membrane (small ao) which is rather pulled 
down as a whole instead of being locally deformed. The 
membrane mediated interactions between bonds then de- 
cay on the length scale ^, and an effective description in 
terms of the bonds alone becomes feasible. 



III. SQUARE LATTICE: NEAREST-NEIGHBOR 
APPROXIMATION 



For a perfectly fiat membrane with ctq = (implying 
X = 0) there is no coupling between bonds, i.e., the bonds 




FIG. 2: Decay of the effective membrane-induced interactions 
m{r) between bonds from Eq. (141. The dotted vertical line 



marks the bond distance a used for the nearest-neighbor ap- 
proximation. 



are independent. The partition sum is then calculated 
easily as 



N 



i=l hi 



1 + & 



/3(eb-et) 



N 



The mean bond density 

1 dlnZ 



o/3(eb-et) 



(16) 



(17) 



shows a continuous crossover between a bound (for eb ^ 
et) and an unbound (for eb ^ et) state with (0) = 1/2 
for Cb = et. 

For a more quantitative analysis of a fluctuating mem- 
brane with (To > we exploit the fact that the effective 
interactions decay exponentially fast, see Fig. [2] There- 
fore, we can approximate the bond interactions by tak- 
ing into account only nearest-neighbors. Moreover, we 
assume the ligand-receptor pairs to be arranged on a 
square lattice with spacing a. For an upper bound of 
the validity of this regime, we choose a lattice spacing 
a = 1^ with nil = m(a) ~ 0.42 and interaction energy 
v ~ etX"^i- Exploiting the analogy with the lattice gas 
we can immediately specify the equilibrium phase dia- 
gram for independent energies eb and et parametrized by 
X, see Fig. |3] Using that at phase coexistence the chem- 
ical potential is /Ltco = — 4i/ and that the critical points 
obey Pv* ~ 2/2.269 [55], the line of critical points is 
determined as 



mi 



-(l-f4mi) 



(18) 



with slope one, where the offset depends on mi. Spec- 
ifying X selects a single critical point. At every critical 
point, a coexistence line eb(et) = zet with slope 



z = 1 — x(l + 4mi) 



(19) 
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FIG. 3: Set of equilibrium phase diagrams for bonds on a 
square lattice with nearest-neighbor interactions (mi ~ 0.42). 
Critical points (filled circles) fall on a line (thick line). The 
shaded region indicates the accessible range of coexistence 
(solid) and crossover (dashed) lines, and the relevant specific 
binding energies. In Sec. |IV| we investigate the dynamics of 
two state points marked by the open circles for x ~ 0.1. 



for et > el terminates. The straight continuation of the 
coexistence Hne below the critical point et < el marks 
the continuous crossover between bound and unbound 
membrane where {4>) = 1/2. 

We restrict our discussion to typical values for the spe- 
cific binding energy eb ranging from bk^T to SSfceT [HJ 
[29] as indicated in Fig. [sj The critical point correspond- 
ing to the upper limit implies Xi — 0.05. Hence, in the 
range ^ x < Xi for eb ^ SS/cbT' the qualitative pic- 
ture for X = persists and a continuous crossover is 
observable. For x > Xi the critical point enters the rele- 
vant parameter region and the first order transition be- 
comes accessible. Increasing x further, the critical point 
moves down the line given by Eq. ( 18 1 to smaller spe- 



cific binding energies, eventually reaching the lower limit 
ej^ ~ SfceT for x — 0.2. Hence, the shaded area in Fig. [3] 
marks the set of coexistence and crossover lines of the 
effective lattice gas that is compatible with our initial 
assumption of a small cooperativity parameter x- A gen- 
uine coexistence between closed bonds and open bonds 
requires both non-negligible membrane fluctuations and 
rather large specific binding energies eb and tether en- 
ergies et. For either weak fluctuations or weak binding 
energies only a continuous crossover between these two 
states is obtained. 

Pushing the analogy with the lattice gas further, in the 
first order transition regime we estimate the free energy 
barrier an initially unbound membrane needs to over- 
come in order to bind to the substrate. From the parti- 



tion sum ( 12 ) we extract the free energy 



N 



(n) 



where the sum in the first term runs over all bond pairs. 
It counts the number of 'broken' links (two neighbor- 
ing bonds in different states) adding up to the interface 
length between bound and unbound domains. Assum- 
ing a circular domain of size n in the spirit of Becker- 
Doring [30], the change of free energy associated with 
this domain is 

AG(n) ~ 2v\/Tm — (eb — zet)n. 

The critical size of the nucleus and the barrier height are 

„t AG* ~ , 



(eb - zetY 



Cb - zet 



respectively. The nucleation time then grows exponen- 
tially with /3AGi. 



IV. DYNAMICS 

Going beyond static properties, in this section we com- 
pare the effective bond lattice gas with the full model for 
the adhesion dynamics close to coexistence. Integrat- 
ing out the membrane degrees of freedom, one expects 
the description of bond dynamics through a Markovian 
stochastic process without memory to be a good approx- 
imation in the limit where the formation of bonds is slow 
compared to the time scale of membrane undulations. 



A. Binding and unbinding rates 

We start with the implementation of dynamics in the 
full model. In Fourier space, the membrane modes obey 
the Langevin equation 



dth, 



't"-q 



-A 



dU 



+ c 



q' 



(20) 



where Aq are the Onsager coefficients that take into ac- 
count the fluid surrounding the membrane [13j and Cq is 
Gaussian white noise obeying the fluctuation-dissipation 
theorem. Further details of how the dynamics of the 
membrane fluctuations is implemented in the full simu- 
lation scheme are discussed in Ref . [HI [32] . 

After a membrane step the heights h{ri) are extracted 
via discrete Fourier transformation and the bonds are up- 
dated. The membrane exerts a force fi on a closed bond 
that is given by the effective spring constant k times the 
separation. The unbinding rates for a bond to rupture 
are then modeled as the usual Bell rates |33|, 



uoff 



Woe 



Pf"'", U = k[h{v,)-h] 



(21) 



with bare dissociation rate and effective size of the 
binding potential set to o-q- The binding rates follow 
through the detailed balance condition as. 



kfe 



■i/3fe[/i(ri)-io]^-|-/3eb 



(22) 
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FIG. 4: Membrane height profile and bonds for a single run of an initially unbound membrane: when the critical nucleus has 
formed (left), when the bond density has reached ~ 1/2 (center), and in equilibrium (right). A closed bond is drawn light, 
an open bond is dark. Parameters as in Fig. Isb). 



These rates are, therefore, independent of the state of 
other bonds. Cooperative dynamics between bonds, e.g., 
the assistance of bond formation, are completely due to 
interactions mediated by the membrane. In Fig. [4] a 
representative example of the membrane nuclcation dy- 
namics of the full model is shown. 

For constructing the dynamics in the effective lattice 
gas, we do not have the explicit height information h{ri) 
of the membrane anymore. The most naive approxi- 
mation would be to assume a constant unbinding rate. 
Somewhat more realistically, we can construct an ap- 
proximation of the membrane force fi = kli({bi}) based 
on the configuration of neighboring bonds. As a first 
implementation we will employ a simple linear interpo- 
lation li = of the bond separation depending only 
on the bond density at a given time. The typical 
height of the unbound membrane is hg. Force balance 
between tether force and non-specific interactions for a 
single closed bond implies 

a'^-/{hi - ho) + kihi - Iq) = 

for the typical height hi of the bound membrane. Linear 
interpolation then leads to the expression 



;((/)) ho{l - 0) + hi<j) -Iq = [ho - lo) 



1 - 



k(j) ' 

0^7 + k 



for the bond separation. Hence, for the unbinding rate 
we obtain 



r,LG,off 



(0) 



Wo ° exp ■ 



nil 



1 - 



32x 



32x 



(23) 



in the reduced quantities. The rate determines the 
time scale of the lattice gas. Since the membrane undu- 
lations slow down the bond formation in the full model, 
this rate will be smaller than wq. The binding rate again 
follows through detailed balance 

^LG,on _ jJ.G,oS^-puY.j(i)bj-Pn ^24) 

and thus depends on the state of the neighboring bonds. 



B. Numerics 

For the numerical simulations we set the bending 
rigidity /3k = 80 and the non-specific strength /37 = 
10^^ nm^^. The bond distance is a = 80 nm with 
bond strength pk = 2.25 x 10~^nm~^. These values 
lead to height fluctuations ctq — 2.1 nm, decay length 
^ ~ 53.3 nm, and slope z ~ 0.73 of the coexistence line. 
The cooperativity parameter is x — O-l- 

In Fig. [5] binding dynamics in the nucleation regime 
of an initially unbound membrane are compared to the 
bond lattice gas employing kinetic Monte Carlo moves 
for two different state points (indicated in Fig. |3]). For 
Fig. [5^) and b), we choose /?et — 28.4 (corresponding to 
I3y = 1.2, ho-lo ~ 50.3 nm) and /3eb = z/3et-h0.2 ~ 21.0. 
For Fig. [HJ:), we move closer towards the critical point 
with /3et ~ 23.7 {I3v ^ 1, ha ~ Iq ~ 45.9 nm) and 
/3eb = zf3et + 0.27 ~ 17.6. The plotted mean densities 
(p{t) are averaged over multiple runs, where single runs 
are shifted such that at t = the critical nucleus has 
formed. Fig. [5|d) shows a system of 256x256 receptor- 
ligand pairs compared to a 128x128 grid in Fig. [5^). 
The curves are qualitatively the same for both system 
sizes, however, the waiting time (not shown) for the crit- 
ical nucleus to form is much smaller in the bigger system, 
which is to be expected. The fluctuations of the bound 
membrane are larger in the full model compared to the 
lattice gas. For comparison, in Fig. [5]d) we show the lat- 
tice gas dynamics for constant unbinding rate fc'"'^'°^(0). 
While it follows the curve of the full model for (j> < 0.5, 
it misses to capture the slowing down of the dynamics 
for higher bond density. The mean-field rates k^'^'°^{(j)) 
perform much better but still for </) > 0.85 the lattice 
gas dynamics for the almost bound membrane becomes 
faster compared to the full model. Moving towards the 
critical point, dynamics of bond formation slows down as 
is clearly visible from Fig. [5]:). 

The bare dissociation rate in the full model is set to 
Wo = 5 X 10^ s^^. This rate is close to the inverse re- 
laxation time of membrane undulations on the length 
scale of the bond distance for the unbound membrane, 
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FIG. 5: Dynamics of the lattice gas (LG) compared to the full model for x — 0-1 (parameters specified in the main text). Shown 
is the density (j){t) for an initially unbound membrane in the nucleation regime close to the coexistence line with fiet ~ 28.4, 
Pexy ~ 21.0 on a) 128x128 and b) a 256x256 grid, c) Closer to the critical point with /3et ^ 23.7, /3eb ^ 17.6 on a 256x256 
grid. The arrows mark the nucleation barrier <j>^ = n^/N. 



Tun ^ rjo? I K ^ 10^^ s with the viscosity of water 77 ~ 
10~^ Pa • s. Such a correspondence imphes that we test a 
dynamic regime where there is no large time scale separa- 
tion between membrane undulations and bond formation. 
Nevertheless, the binding dynamics in Fig. [5] show a very 
good agreement between full model and lattice gas. Re- 
alistic dissociation rates are typically somewhat smaller 
and, therefore, in a regime where the lattice gas will per- 
form even better. The rates have been extracted 
from scaling the time evolution to fit the full model 
and have been determined as: Wq'^ /wq ~ 1.8 x 10"'^, 
2.8 X 10-3, 2.1 x 10"'' for Fig. [5]a)-c) , respectively. 



V. CONCLUSIONS 



costly simulations of the full model with the fast lattice 
gas simulations reveals a very good qualitative agreement 
in the increase of the average bond number as a function 
of time. We, therefore, expect our model to be useful 
in more complex studies of, e.g., moving ligands and the 
influence of defects on membrane adhesion. 
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Beginning with a Hamiltonian that incorporates the 
essentials of specific adhesion to a flat substrate, i.e., 
bending energy of the membrane, a non-specific adhesion 
potential, and the energy contribution of receptor-ligand 
pairs, we integrate out the membrane shape in the par- 
tition function of the system. In the limit of small mem- 
brane roughness and/or weak stiffnesses of the tethers 
attaching the receptors to the substrate, our approach 
leads to an effective free energy that is isomorphic to 
that of a lattice gas model used to describe liquid-gas 
transitions. We are able to give an explicit expression for 
the effective interaction potential between bonds. When 
the bond sites are far enough apart it is sufficient to con- 
sider only nearest-neighbor interactions. For bonds on 
a square lattice the equilibrium phase behavior follows 
immediately. The analogy to lattice gas models also ex- 
tends to nucleation theory, yielding expressions for the 
critical nucleus size of an adhered membrane patch and 
the energy barrier that needs to be crossed. 

We further show that the lattice gas model may also be 
utilized to study the dynamics of the adhesion process. 
To this end, we compare simulations of the full Hamil- 
tonian with a kinetic Monte Carlo scheme for the lattice 
gas. The comparison of results of the computationally 



Appendix A: Membrane fluctuations 

We expand the height profile /i(r) into Fourier modes, 

/^(r)=^V-'i-, (Al) 
q 

where /i* = /i_q. For a square system of size L, the ac- 
cessible wave vectors are q — (2tt / L){nx, riy) with integer 
rix, TLy, and q = |q|. The Helfrich energy M becomes 

In the non-specific interactions ^ 

we split off the zero-mode h = ho contribution. 

We want to integrate out the fluctuations of the mem- 
brane under the constraint that at the height is 
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h + (Tohi. This constraint is expressed as 

N ^ 

X exp |i^ Ai[/i(ri) -h- a^hi]^ 



In terms of Fourier modes we have 

h{Yi) = + ^ ' [/iqe-"J""' + /l- 

q#0 



where the sum runs over independent modes q (for 
a q exclude — q from the sum). SpHtting the coefficient 
/iq = /iq + i/iq iuto real and imaginary part, the integral 
over all modes (except the zero-mode) becomes a product 
of two independent integrals. These are simple Gaussians 
and read 

Z'=( Y[Ah'^ 

X exp J - ^ '/3A(^g4 + ^)(K^f + i ^ 

for the integration over real parts. Analogously, Z" is 
obtained through replacing h!^ by /iq and Cq by Cq. The 
coefficients are 



N 



N 



Cq = 2^ AjC0s(q • r^), Cq = 2^ AjSin(q • r^). 



i=l 

The result is 



j=i 



Z'Z" exp I -io-g ^ TOij AiAj I 



with coupling matrix 



E 

q#0 



,cosq ■ (r, - rj) 

Kq* + 7 



We notationally suppress the non-exponential prefactors. 
The integral over the A^ implementing the constraints 
again is a simple multi-dimensional Gaussian integral, 



I 



(27r) 



N 



dAi---dA 



N 



X exp I -^ctq E ^^o^i^j ^ i E '^o/ijA, | 
[(27r)^ dct m]-^'^ exp ^(to-1),j/i,/Ij | 



The remaining energy is a function of the mean height, 



E(h) =^Ch~ Sho)' +Y^bJ'lCh + aoh,f - Ebj , 

i=l 



where we have shifted /i — > h + Iq for convenience and 
defined SHq = Hq — Iq. We finally integrate over the mean 
height 



with 



1 (fcg^o E»^»^» -7A(5feo)^ 

2 7^4 fcAf0 



7A 



Appendix B: Effective interaction 



To obtain a more useful expression for the bond inter- 
actions, we approximate the sum over modes in Eq. ([t]) 
by a two-dimensional integral, 



m{r) = 



2^^PKal Jo Jo 1 + m 



I q cosfqr cos 09] 



27r/3Kcr2 



dq 



"'0 



qJpjqr) 



where ^ = (k/7)^^^ and Jo(a:) is the zero-order Bessel 
function of the first kind. Performing the integral [34] 
leads to the effective interaction 



4 

= -^nn — 7^^2 keio(?-/0 = keio(r/^). 

27r/3K(cro/C) TT 



Since keio(O) = — 7r/4 we have to(0) = 1. Here, keio(a:;) is 
a Kelvin function defined as 



keio(x) EE luiKoixe'^''/^), 



(Bl) 



where Ko{z) is the zero-order modified Bessel function of 
the second kind. 
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